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Abstract 

The structures of confining vortices which underlie pure SU{3) Yang-Mills theory 
are studied by means of lattice gauge theory. Vortices and monopoles are defined 
as dynamical degrees of freedom of the gauge theory which emerges by center 
gauge fixing and by subsequent center projection. It is observed for the first time 
for the case of SU{3) that these degrees of freedom are sensible in the continuum 
limit: the planar vortex density and the monopole density properly scales with the 
lattice spacing. By contrast to earlier findings concerning the gauge group SU{2), 
the effective vortex theory only reproduces 62% of the full string tension. On the 
other hand, however, the removal of the vortices from the lattice configurations yields 
ensembles with vanishing string tension. SU{3) vortex matter which originates from 
Laplacian center gauge fixing is also discussed. Although these vortices recover the 
full string tension, they lack a direct interpretation as physical degrees of freedom in 
the continuum limit. 
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1 Introduction 



Since QCD was recognized as the theory of strong interactions by means of high energy 
scattering experiments, the question arose as to whether QCD also explains the absence of 
quarks in the particle spectrum. After extensive numerical simulations had become feasible 
with modern computers, it became clear that the pure gluonic theory already bears witness 
to quark confinement: the static quark anti-quark potential rises linearly for large distances 
due to the formation of a color-electric flux tube [Tj. Moreover, the low-energy model, in 
which a fluctuating bosonic string plays the role of the effective degree of freedom, predicts 
a characteristic 1/r correction to the potential at large distances. This picture recently 
received viable support from lattice simulations which verified the dependence of this term 
on the number of dimensions at a quantitative level 0. A major challenge of modern 
quantum field theory is the question: Why does the color-electric flux tube form? 

Over the recent past, lattice gauge simulations have strengthened the idea that topological 
degrees of freedom, which are characteristic for the non-Abelian nature, are relevant for 
confinement. Among those, color- magnetic monopoles and center vortices are under intense 
discussion (for a most recent review see Here, we will focus on the vortex picture of 
confinement. 

The central idea is to simplify pure Yang-Mills theory under the retention of its confining 
capability, hoping to filter out degrees of freedom which meet two criteria: (i) the degrees 
of freedom are sensible in the continuum limit of lattice gauge theory, and (ii) they are 
closely related to confinement of pure Yang-Mills theory. 

Gauge fixing and projection techniques have proven to be convenient for these purposes. 
Center gauges have been designed to maximize the importance of center degrees of freedom, 
and the projection SU{N) — > was proposed for the simplification process [U 
Vortices and monopoles appear as the dynamical degrees of freedom of the Zn gauge 
theory. Indeed, criteria (i) and (ii) were seen to be satisfied in the case of a SU (2) gauge 
theory if the maximal center gauge (MCG) is adopted jH El El El • 

Of particular interest is the case of the SU{3) gauge group because of its relevance to 
the theory of strong interactions. In the present paper, a thorough study of the confining 
vortices is performed for the gauge group SU{3). For the first time the planar vortex 
density as well as the density of the Z^ monopoles are reported to properly extrapolate 
to the continuum limit in the case of MCG. In sharp contrast to the case of SU{2), the 
vortices fail to recover the string tension to its full extent. We will see below, however, that 
there is still a close relation of the MCG vortices to confinement: removing the vortices 
yields a model theory with vanishing string tension. These findings will be contrasted to 
those obtained in the Laplacian center gauge. 

The techniques for extracting static quark potentials as well as the numerical setup are 
explained in the next section. Details of the definition and construction of the vortex 
matter are given in section |21 To which extent the vortex matter is able to reproduce the 
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SU{3) string tension is studied in section^ Thereby, new high-precision data for the case 
of SU{2) are presented. These data serve as a "contrast agent" for the findings for the 
SU{3) case. In sectional the properties of vortices and Z3 monopoles are discussed in the 
continuum hmit. Conclusions are left to the final section. 

2 The static quark potential 

Focal points are the study of the relevance of the vortices for the SU (3) string tension and 
of the properties of the vortex matter emerging in the continuum limit. For these purposes, 
a careful determination of the physical value of the lattice spacing is mandatory. This is 
the subject of the present section. 

2.1 Numerical Setting 

Results of the simulation of pure SU{N) gauge lattice gauge theory will be presented 
below for N = 2 and N = 3, respectively. The dynamical degrees of freedom are the 
unitary matrices Ufj^{x) e SU{N). Configurations {U^{x)} will be generated according to 
the Wilson action 

^ = ^ E ^tr(p,.(a;) + Pj,(a;)) , (1) 

P^.u{x) = U^{x)U,{x + fi)Ul{x + u)Ul{x), (2) 

where P^v{x) is the plaquette. The update is performed using the Creutz heat-bath algo- 
rithm |H5 for the case of SU{2). Updating the diagonal SU{2) subgroups as proposed by 
Cabibbo and Marinari is performed in the case of SU{3). Each ten heat-bath sweep is 
accompanied by four micro-canonical reflections in order to reduce autocorrelations. Mea- 
surements were taken after twenty such blocks of sweeps. Most of the data are taken on a 
L'^, L = 16 lattice. In the case of the SU{2) gauge group, lattices with L = 24 were also 
studied. 

The static quark anti-quark potential V{r) will be extracted from planar, rectangular 
Wilson loops W{R,T) of extension Rx T, i.e., 

W{R,T)^ oc exp{-V{r) aT} , r := Ra, T : large, (3) 

where a is the lattice spacing. 

2.2 Overlap enhancement 

In order to extract the physical signal from noisy Wilson loops, so-called overlap enhance- 
ment has been proven to be an important tool fOl m For the case of S'f/(2)^, we 

^For the case of SU{3), a novel method will be proposed below. 
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Figure 1: Linearity of the function — In {W{R, T)) in T for several values of R (left panel). 
The static quark anti-quark potential for SU{3), 16^ lattice and /3 = 5.9 (right panel). 

closely follow the procedure in and define the cooled spatial links by 

Uiix) = VnY. ^'^(^) ^^(^ + ^) + ^) ' ^ = 1, 2, 3 , (4) 

k 

where k runs from —3 ... 3 and k = i is excluded from the sum. Vn is the projector onto 
the 'closest' SU{N) element. In the case of SU{2) and for M = oq + ira, being the 
Pauli matrices, the effect of the operator V2 is 

V2M = -^^= (ao + ^rd) . (5) 
+ ^ ^ 



Details concerning V3 can be found in [12]. Temporal links are unchanged, i.e., 

no(x) = Uoix) . (6) 

Since the projection Vn, = 3 is "expensive" from a numerical point of view, we used 
a different method to define the cooled spatial links Ili{x), i = 1,2,3. Let us define the 
action of the spatial links of a given time slice t by 

j>fc,1...3 

where Pik{x) ^ is the plaquette calculated from the spatial links. In addition, we define 
three different embeddings of the SU{2) matrix aQ + ira, Og + a^ = 1 into the group S'f/(3), 
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Figure 2: The static quark anti-quark potential for the SU{2) gauge group (left panel) 
and for the case of the SU{3) gauge group (right panel). 



I.e., 



'10 

V^^^ = I ao + as 02 — iai 
,0 — a2 — iai ao — 



00 + 03 a2 — iai 
V^^^ = ( 10 

-a2 — iai ao — 03 



ao + as a2 — iai 
V^^"^ = I — a2 — iai c^o ^ % 
1, 



(9) 



Let us now consider a particular spatial link Ui{x). Substituting U' := V^^^Ui{x), we locally 
maximize the action S'(3) with respect to V^*^^-*. Subsequently, we replace U' by U" = V^'^^U' 
and maximize with respect to V^'^\ and setting U'" = V^^^U", V^^^ is chosen to maximize 
S'(3). Finally, we define 

Uiix) = V^^^ y(l) jj^^^-^ _ (10) 

We then visit the next link on the lattice. One sweep has been performed when all spatial 
links of the lattice have been visited. The advantage of the present procedure is that the 
maximization of ^(s) with respect to one of the SU (2) subgroups can be implemented very 
efficiently. 

In order to achieve a good overlap with the groundstate, the above procedure for determin- 
ing the links n;(a;) is applied recursively, and the Wilson loop expectation value {W{R, T)) 
is calculated from the configurations {n^(x)} rather than the ensembles {U^{x)} . The 
average of ^(s) over the time slices divided by the total number of spatial links serves as 
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0.26(2) 


0.24(1) 


0.146(3) 


0.145(2) 


0.0752(7) 


0.0754(5) 


0.0391(3) 


0.0373(1) 



Table 1: String tension a in units of the lattice spacing a for the case of the SU{2) gauge 
group. 

litmus paper for the overlap. It turned out that ten sweeps are enough to yield more than 

0. 99 ground state overlap. Good overlap is also signaled by the quantity 

-\n(wiR,T)y (11) 

which already shows a linear behavior in T for T > 3. This is illustrated for the case of 
a SU{3) gauge group, L"^ = 16^, and P = 5.9. The final results are obtained from 100 
independent measurements. The symbols in figure [U left panel, represent the lattice data 
of the quantity (|lip: the lines are linear fits in T, i.e., 

-In (W{R,T)^ = + 6 . (12) 

The coefficients 7 can be interpreted as V{r)a. The latter quantity is shown in figure ^ 
right panel (symbols). The line is a fit according to the function 

V{r) a = cR - ^ + Vo , r = Ra , (13) 
R 

where the parameter c can interpreted as the string tension in units of the lattice spacing, 

1. e., c = cra^. This method for the calculation of the string tension was used, e.g., in [T^ . 
For the present example, we find 

aa^{p = 5.9) = 0.0701 ± 0.0005 , (14) 

which is in agreement with the value reported in the literature, i.e., a (/5 = 5.9) = 0.073 

[HI eg. 

2.3 The scaling relation 

In order to mark a dimensionful quantity as a physically sensible one in the continuum 
limit a — > 0, it is crucial to express this quantity in units of a physical reference scale. 
Throughout this paper, the string tension a will serve this purpose. 

It is time consuming but straightforward to determine the dependence of a a^(/3) on j3. 
The results for the gauge group SU{2) are summarized in tabled Those for the case of a 
SU (3) gauge group are shown in table |21 

These values are in good agreement with those reported, e.g., in [T^ . 
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0.32(1) 


0.26(1) 


0.169(3) 


0.104(1) 


0.0701(5) 


0.0514(3) 



Table 2: String tension a in units of the lattice spacing a for the case of the SU{?>) gauge 
group. 

3 The vortex texture 

Concerning the criteria (i) and (ii) of the introductory section in^ one firstly notes that 
the properties of the SU{2) vortices extrapolate properly to the continuum limit 6 if the 
maximal center gauge (MCG) is used. Moreover, the MCG vortex theory correctly repro- 
duces the de- confinement temperature The phase transition acquires a geometrical 
picture: is appears as a vortex de-percolation phase transition fUlEj which already points 
towards a weak vortex interaction (see also JH])- Secondly, it was observed that this ge- 
ometrical picture correctly reproduces the finite size scaling of the 3D Ising universality 
class pn]- This shows that at least for temperatures close to de-confinement, MCG vortices 
interact weakly. 

It turns out that the non-locality induced by gauge fixing is crucial for property (i) above: 
it was analytically shown that if unfixed lattice configurations are projected onto vortex 
configurations, the complete static quark potential (including the Coulomb term) is ob- 
tained |2ni- At the same time, the properties of this vortex matter strongly depend on 
the size of the lattice spacing. Thus, these vortices lack an interpretation in the contin- 
uum limit. In these respects the non-locality of the approach to the vortex matter is an 
advantage. On the other hand, however, this non-locality generically makes it difficult 
to access the vortex matter in practical simulations: using a variational gauge condition 
as, e.g., in the case of MCG [H], the definition of the vortices is ambiguous as a result of 
the inability to localize the global maximum of a non-linear functional. At least for small 
lattice sizes, the so-called Gribov ambiguity might have a significant influence on physical 
observables fTHU^. 

Marked progress concerning the Gribov problem was made with the construction of the 
Laplacian gauges ^j: the non- locality of the gauge fixing is preserved while the maximiza- 
tion is replaced by an eigenvector problem, which can be handled by present-day algorithms. 
The Laplacian version of the MCG was firstly proposed for a SU{2) gauge group in [21] 
and generalized to the SU{3) case in [25j. A further improvement was reported in pU] . 
Vortex matter of the Laplacian center gauge (LCG) is unambiguously defined and recovers 
the asymptotic string tension for both gauge groups, SU{2) and SU{3). In the case of 
a SU{2) gauge theory, it was, however, observed that the LCG vortices are produced in 
large abundance, implying that they lie dense (nevertheless in a controlled way) in the 
continuum limit [7j. This indicates a rather strong interaction of LCG vortices, which 
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might render it difficult to mimic LCG vortex matter in a low-energy effective model. 

The present section brieffy reviews details of the center gauge fixing procedures and the 
vortex projection techniques with an emphasis on the case of SU{3). 

3.1 The ideal center vortex cluster 

In order to reveal degrees of freedom which are relevant for confinement, we are looking 
for configurations {Zfj_{x)}, Z^{x) G Zn which best represent the full link configurations 
{f/^(x)}, U^{x) G SU{N). Thereby, Z^ represents the center of the group SU{N), i.e., 

r 27r ^ N N 

Z^{x) = expji— m| , - Y < m < y , (15) 

where m is an integer. There is an optimal choice of the gauge for which the overlap of the 
center configurations with the full ones is maximal. Let us denote the gauge transformed 
links by 

u^{x) = n{x) u^{x) n^{x + fi) , n{x)ESU{N). (le) 

In oder to obtain the ideal center configurations, we minimize the functional 



it 



Ul\x) - Z^{x) min (17) 



with respect to fl{x) and Z^{x). The selection of fl{x) implies the choice of a gauge. We 
will call this gauge ideal center gauge (ICG) throughout this paper. The condition (fTTj) 
directly implies that the overlap, i.e., 

is maximized. Ni is the number of links of the lattice, and — l<i?<l. R = 1 means that 
the link configuration {Ufj,{x)} can be entirely expressed in terms of center elements after 
a suitable gauge has been chosen. 

The maximization of R with respect to the center elements Z^{x) can be performed locally: 
with 

1 27r 
— trf/^(x) = ui exp{V}, Z^{x) = exp{iip^^} , cp^ = — m , (19) 

where / = {x, fi} specifies the link, one finds 

The optimal choice is obtained by choosing m (fT^ in such a way that is closest to ipK 
This mapping 

SU{N)^Z^: U;!{x)^Z,{x) = Z(u^{x)) (20) 

is called center projection. Inserting Z^{x) from (j20|) into (jl8|) . the overlap R must then 
be maximized with respect to the gauge transformations, i.e., Q{x). 
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3.2 Center gauges 

Since the mapping (j20|l by no means depends smoothly on Q{x), an iteration over-relaxation 
algorithm which iteratively determines Z^{x) and VL{x) from ()18|1 can hardly work. State 
of the art would be to determine the desired quantities Zfj_{x), VL{x) by the technique of 
"simulated annealing". However in this case, determining, e.g., Vl[x) to the precision which 
is needed for the vortex analysis is extremely costly from a numerical point of view. This 
approach is beyond the scope of the present paper and is left to future investigations. 

In order to make the gauge fixing efficient by means of iteration over-relaxation, we assume 
that in the optimal case, iiU^{x) comes close to a center element. In this case, we 
relax the condition which constrains Z(U^{x)) in (I^U]) to a center element. There are two 
possibihties for doing this: 



z[u^;{x] 
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oc 



oc 



iiU'Ux)^ 



Hence, we find the gauge conditions 
R 
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(21) 
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X,fl 



both have been advertised in the literature ^7\. It is not clear at the beginning which 
of the above possibilities yields the larger overlap (jl8p . It might turn out that even a 
background-field-dependent admixture of both possibilities (j21l22j) is best for the present 
purposes. Since the algorithm for the so-called "mesonic" gauge was already studied 
in the literature to a large extent j2Zl I2H| , we will employ the gauge condition Rmes — ^ max 
for the determination of Q{x) and perform subsequent center projection along the lines 
outlined in the previous subsection. 



Even once we have agreed on one of the suboptimal gauge conditions ()23I24|) . there is still 
the problem of the Gribov ambiguities: since the overlap R is a non-linear functional on 
Q{x), detecting the global maximum of R is practically impossible for reasonable lattice 
sizes. The choice of the local maximum of R which is implicitly defined by the algorithm 
determines the gauge. Although the physics extracted in the latter gauge can be highly rel- 
evant for confinement, the definition of gauge by the numerical procedure is unsatisfactory. 
The Gribov problem can be solved by adopting the Laplacian version [2^1211 1231^01 of the 
maximal center gauges. A detailed study of these gauges and the corresponding vortex 
matter can be found for the case of SU{2) in [7] and for the case of SU{3) in [22]. Here, we 
will briefiy outline the procedure for the case of SU{3) and refer the reader to reference j2Sl 
for details. 
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The generators t°- of the SU{N) algebra satisfy the equation 

C = Smk - ^S^m S^k) , tl t^t' = ^S'^' . (25) 

With the help of this identity, the "mesonic" functional Rmes f!23|) can be written as 

R^es = ^ E [tr 0^(x) R^{x) 0{x + fi) + l] , (26) 

X,fl 

where the adjoint matrices are defined by 

0''\x) = 2tr|r n{x) n\x)^ , Rf{x) = 2tr|r U^{x) [/^(x)} . (27) 
Equation ()26|) can be written as 

Rmes = ^Tr 7^ O + 1 , (28) 

where we have introduced, e.g., the vector O of the combined coordinate and color space, 
{0°'^{x)} — > O. Up to a term proportional to the unit matrix, TZ is the adjoint Laplacian 
operator, i.e., 

K = I ^V'^-^f^ + ^?(^ - '^^.-M 1 • (29) 



Note that the vector O is subjected to the constraints that the set of vectors with 
Oab{x) = {n^{x), ri^^x), n^{x)}ab is orthonormal. At the heart of Laplacian center gauge 
fixing, one relaxes these constraints and seeks the — 1 largest eigenvalues of the super- 
matrix TZ. These tasks can be unambiguously performed with present-day algorithmic 
tools. From the corresponding eigenvectors, the adjoint gauge transformations 0"*(x) 
are reconstructed at each site with the help of Gram-Schmidt orthogonalization. Abelian 
monopoles and vortices appear as defects in the latter step of reconstructing the gauge 
transformation. Technical details of this gauge fixing are presented in j2^. Finally, we 
point out that the Laplacian gauge also seeks to maximize the "mesonic" gauge condition. 
However, the re-orthogonalization of the eigenvectors implies that the values for the overlap 
R (jlSp are significantly smaller than that achieved by maximizing Rmes ()23p with the help 
of an iteration over-relaxation procedure. 

3.3 Identifying vortex matter 

In order to reveal the vortex matter of Zn gauge theory, we define 

v{p) ■= n^/) l = {x,f^}, v{p) e Zn , (30) 

l€p 
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where p = {x, ^ < u) defines an elementary plaquette on the lattice. One says that a 
vortex of center charge z pierces the plaquette p if 

v{p) = z , z = expji— ^1 , (31) 

where — y < Lp{p) = ^^„{x) < y is called the center flux. For SU{N > 3), one defines the 
conserved Z]\f monopole current by 

mf,{x) = ^ ef,„ap (faisix) , A^m^(x) = 0, (32) 

where 

A^^{x) := $(x + /i) - ^{x) . 

In order to reveal the span of the monopole charge m^, we consider an elementary, spatial 
hypercube c = (x, a < /3 < 7). One easily verifies, using the Abelian nature of the group 
Zjv, that 

1 = ]^£'|^t'(p) = expjz^ eifcmAj(^fcm(a;) I , (33) 

where the sum over i, k, m runs from 1 ... 3, and D[v{p)] = w^(p) if the normal vector of 
the plaquette is anti-aligned with the normal vector of the relevant surface element of c, 
and D[v{p)] = v{p) otherwise. Hence, the monopole charge comes in integers, i.e., 

7714 = ^ eikmAiipkm{x) = k , k : integer . (34) 

Finally, we point out that there are no center monopoles in the case of a SU{2) gauge 
group, i.e., m^{x) = 0. 

It is convenient for model building to define the vortex matter on the dual lattice, where 
the link I, the plaquette p and the cube c is mapped onto 

/ — > c* , p — ^ p* , c — > r . (35) 

The vortex field of the dual lattice is defined via the identification 

v{p*) = v{p) , for p — > p* . (36) 

The identity ()33p can be transformed into an identity for dual fields only: 

1 = n ^(p*) • (37) 

The latter equation implies that the vortices either from closed world sheets on the dual 
lattice or, for SU{N > 3) only, multiples of vortex world sheets merge at a closed 
monopole trajectory. 
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Figure 3: The static quark anti-quark potential for the SU{2) gauge group: full ensembles, 
vortex projected ensembles, and ensembles where the vortices have been removed (see (jHSI))- 
Lattice size: 24^. 

4 Dominance of the static quark potential 

Preliminary evidence that the MCG vortex matter recovers the string tension of pure 
SU{?>) gauge theory was presented in [2Zj. The so-called indirect center gauge, i.e., max- 
imal Abelian gauge fixing with a subsequent fixing of the center gauge, was investigated 
in lin]. There it was observed that the string tension obtained from Abelian monopole 
configurations as well as from Z3 vortex ensembles is significantly smaller than the string 
tension of complete SU{'i) gauge theory. On the other hand, it is interesting to note that 
the SU{?>) string tension is recovered to full extent in the case of the LOG |25]. 

Here, we will see that the MCG vortices act like the vortices defined in the indirect center 
gauge: roughly 62% of the full string tension is found. We will therefore find that the 
MCG vortex ensembles genuinely differ from those of the LCG. 

4.1 The case of a SU{2) gauge group revisited 

In order to contrast the findings concerning the SU (3) gauge group, which will be presented 
below, with the findings for the case of SU{2), we briefiy discuss our numerical results for 
the latter case in this subsection. 

The maximal center gauge (j23p is implemented by using the iteration over-relaxation pro- 
cedure which is described in detail in This procedure defines the gauge. The corre- 
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75 
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35 


55 



Table 3: Number of independent SU(2) configurations used for the calculation of the 
static potential (see figure E))- 

sponding vortex degrees of freedom are defined by the projection (j2(Jj) . which becomes in 
the present case: 

SU{2) ^ Z2 : U^{x) ^ Z^{x) = sign tr f/^^(x) . (38) 

It turns out that this procedure produces vortex configurations with sensible properties in 
the continuum limit ^ and with a close relation to the physics of confinement [HllllElEl- 
So far, vortex matter with "best" properties in the continuum limit seems to obtained 
with a 'preconditioning' by performing the Laplacian center gauge (see subsection 13. 2|) and 
subsequent maximal center gauge fixing [2ll2ni- This approach also alleviates the Gribov 
problem, but its implementation is numerically "expensive" . 

In order to reveal the relevance of the vortex texture for the physics of confinement, one 
firstly calculates the static quark potential from the vortex configurations. Secondly, one 
defines a toy Yang-Mills theory by 

K{x) = Zlix) U^{x) , (39) 

where the vortex texture has been removed "by hand" from the lattice ensembles. It 
was found [U E] that the vortex configurations reproduce the linear part of the potential 
to a large extent. In addition, the potential evaluated from the modified configurations 
{U'^{x)} has lost its linear rise and shows a Coulomb type of behavior. Both observations 
are summarized by the term 'center dominance of the potential'. 

Figure 01 illustrates these observations using our lattice results. Only data which were 
obtained on a 24"^ lattice are shown. The number of independent configurations employed 
for the calculation of expectation values is listed in table El The full configurations {Ufj_{x)} 
as well as the configurations {U'^{x)} from which the vortices have been removed "by 
hand" are subject to the overlap enhancement described in subsection 12.21 The Z2 vortex 
configurations already possess good overlap with the ground state, and no enhancement is 
used in this case. The line which fits the "vortex only" data in figure El corresponds to a 
string tension of 97.7% of the full string tension, we stress that these findings have been 
obtained with the most naive version of the maximal center gauge (described in 0). 

Finally, we point out that the quality of dominance is affected by the choice of gauge, i.e., 
the Gribov copy [21] , and that the Gribov effect is influenced by the lattice volume [221 ■ 
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Figure 4: The angle <y9*^™^ of the Wilson loop expectation values ()43|) as a function of the 
Wilson loop size n (left panel) and as a function of the size in physical units (right panel). 

4.2 Vortex-limited Wilson loops: SU (3) gauge group 

Let denote a. n x n planar Wilson loop calculated within the particular config- 

uration U^{x). The same object is evaluated with the configurations obtained from 
center projection ()20j) after direct maximal center gauge fixing using the "mesonic" gauge 
condition (|^. The result is called Wn[Z^]. Since Wn[Z^ G Z^,, the latter Wilson loop can 
be characterized by the number m G { — 1, 0, 1}, i.e., 

Wn[Z^] = exp{z|^ m} . (40) 
The expectation value of the Wilson loop, 

Wn = {Wn[U^]) , (41) 

is obtained by averaging over the ensembles {U^{x)}. In addition, we can define expectation 
values 

H^M = (Wn[U,]) , (42) 

where only loops Vrnit^A*] taken into account; the corresponding quantity Wn[Z^] belongs 
to the sector m. Decomposing 

^ exp{z¥,f)}, (43) 

one expects that for large loops 

27]- 

hm = - m , (44) 

n— >oo o 
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if center vortices dominate the Wilson loop expectation value. The latter relation can be 
checked by lattice simulations. The quantity v^rT^ is shown for [3 = 5.6 (227 independent 
measurements) and for (3 = 5.8 (160 independent measurements) in figure 0] as a function 
of n (left panel). It seems that the relation ()44|1 is indeed satisfied for large n. If we plot 
the angle cp^n"^ as a function of the physical size of the Wilson loop, i.e., I = na{(3), we 
observe that the data for /? = 5.6 and for /? = 5.8, respectively, roughly fall on top of the 
same curve (right panel). 

Let us interpret these findings from a random vortex model point of view. Following [2E], 
we assume that center vortex intersection points possess a finite correlation length Lc- 
Thus dividing the minimal area A of the planar Wilson loop into squares of size A > L^, 
the center fluxes ip G {—1, 0, 1} through different squares are essentially uncorrelated. Let 
p(m. A) denote the probability of finding center flux m through the area A; we define the 
"mesoscopic" vortex density by 

Pmes ■■= p{m,A)/A, A = Ti^a^f]). (45) 

Assuming vortex dominance, we might approximate 

Wn = (W^n[f/^]) ~ (w^n[^M]) • (46) 

Using the fact that center fluxes are uncorrelated by construction, one obtains 

^ (Wn[Z,]) = {Z)^/^, (47) 
where (Z) the average flux through the area A is 

(Z) = P(m,^) expjiym}. (48) 



m=—l..l 

Hence, the string tension in the center flux approximation is given by 



A 



5Z exp{zy m} = --^ In (^1 - 3 , (49) 



m=-l..l 



where we have assumed that the center symmetry m -—>■ —m is not spontaneously broken, 
i.e., 

1 

p{l,A) = p{-l,A) =: p{A) , and ^ p{m,A) = 1. (50) 

m=—l 

A particular case is obtained by considering that the vortices which are defined at the level 
of the elementary plaquette are uncorrelated (naive random vortex model). In this case, 
one finds 

p{A = a^) = - <^ 1 , and thus cTc/ ~ - p . (51) 
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Thereby, p is the "microscopic" vortex density, i.e., po? is the probabihty of finding a 
non-trivial center fiux through a given plaquette (no matter whether m = — lorm = l). 

Using the numerical data above, it is possible to estimate the center fiux correlation length. 
Let us define the "half- width" L1/2 by the length of the Wilson loop at which 

¥^^'^^ = ^1/2) = \- (52) 

The findings (see figure El right panel) suggest that 

Lc > Li/2 ^ 0.8 fm, (53) 

where we have used a string tension of cr = (440 MeV)^ as a reference scale. Finally, let us 
check whether the naive random vortex model of uncorrelated vortex plaquettes is realistic. 
For this to be the case, the relation 

f ^ - f Ll;, - ^ L?/, ^ I L?/, « 1 (54) 

must hold. However, one finds, using (j^ . 

I ^ 1-03 , (55) 

implying that the naive random vortex model seems not always to be justified. 

4.3 The "mesonic" center gauge for SU{3) 

In a first step, the "mesonic" gauge condition is installed with the help of the iteration 
over- relaxation algorithm described in detail in [771 1?8|| . The link elements Z^{x) are 
defined by center projection ((201) • the case of SU{2), we will compare the static quark 

potential obtained from full link configurations (see section |21) with the one calculated with 
link ensembles {Z^(x)}. In addition, the toy model is defined by configurations {U'^{x)} 
()39|) from which the vortices have been removed "by hand". From the results of the 
previous subsection, we expect that the string tension is lost in the latter case. Our 
numerical findings using 100 independent measurements are summarized in figure 

We find that the potential calculated from vortex configurations scales towards the con- 
tinuum limit, i.e., the data obtained from different /3 values fall on top of the same curve 
if the V{r) and r are expressed in physical units. In addition, one observes "precocious" 
linearity: the potential is linear even at small distances as is the case for a SU{2) gauge 
group. In contrast to the case of an SU{2) gauge group, the center projected string tension 
is only 62% of the full string tension. The value of string tension (in lattice units) after 
center projection is in agreement with the finding in |27j for a 12^ lattice and (3 = 5.6. In 
the latter article, however, the quoted value of the full string tension is underestimated. 
Using reliable values, the ratio of projected and full string tension is in agreement with the 
findings reported here. 
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Figure 5: The static quark anti-quark potential for the SU{3) gauge group: full ensembles, 
vortex projected ensembles, and ensembles where the vortices have been removed (see (jSHI))- 
Lattice size: 16^. 

On the other hand, removing the vortices (see dSHl)) produces configurations which are 
compatible with a vanishing string tension. There is a subtlety for obtaining this result: 
the lattice volume must be large enough^. It appears that the lattice size of 16^ seems to 
be too small for /3 as big as 6.0. 

Since a removal of the vortices results in a loss of the string tension, even if the vortices 
only amount to 62% of the full one, the question arises as to whether additional degrees 
of freedom which reside in the U{1) x U{1) Abelian subgroup are responsible for the 38% 
string tension completing the vortex contribution. Candidates for such degrees of freedom 
are color magnetic monopoles. To answer this question, we implemented the "mesonic" 
gauge condition and subsequently projected the gauged configurations onto Abelian 
ones: 

U;{x) V,ix) , V,{x) e f/(l) X f/(l) . (56) 
For these purposes, the off-diagonal elements of U^{x) were dropped, i.e., 

U^{^) — U^{^) ■■= diag(f/,^(x)) , (57) 

and V^{x) is given by the SU{3) element which is "closest" to U^{x) (see discussion in 
subsection 12. 2|) . One verifies that indeed V^{x) G f/(l) x U{1). In addition, we investigated 

thank M. Faber for this remark. 
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— » vortex only 
■ U(l)xU(l) removed 
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2 4 6 8 
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Figure 6: The static quark anti-quark potential for the SU (3) gauge group: full ensembles, 
vortex projected ensembles, and ensembles where the vortices have been removed (see (jSHI))- 
Same with U{1) x U{1) projection. Lattice size: 16^, /3 = 5.8. 

ensembles {Ul^^^{x)} which are complementary to the U{1) x f/(l) configurations: 

U^'^'^ix) = Viix) U^ix) . (58) 

Our numerical findings for a 16^ lattice aX (3 = 5.8 are shown in figure El we find that 
the string tension calculated from U{1) x f/(l) configurations is marginally larger than the 
string tension from vortex projected configurations. As expected, configurations U-^^^^x) 
from which the Abelian subgroup was removed do not support confinement. 

These results are compared with those in ^3]: there, versions of the so-called Maximal 
Abelian Gauge were investigated. These gauges are most suitable for a projection of 
configurations U^{x) onto the Abelian subgroup U{1) x U{1). Also in these cases, the 
string tension extracted from U{1) x U{1) configurations is substantially smaller than the 
full string tension. 

4.4 Laplacian center gauge 

The previous subsection showed that the "mesonic gauge" produces vortex matter 
which recovers only 62% of the full string tension. The question is whether Z3 matter 
is able to give the full result for the string tension at all. The answer was already given 
in 1^: vortex matter which is defined by the Laplacian gauge condition (see subsection 
13. 2p reproduces the linear rise of the static quark potential in the continuum limit. Here, 
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Figure 7: The static quark anti-quark potentials for the SU{3) gauge group using the 
Laplacian center gauge: full ensembles, vortex projected ensembles, and ensembles where 
the vortices have been removed (see (jHEj))- Lattice size: 16"^, P = 6.0. 

we briefly report our findings. We investigated the somewhat extreme case of a small 
physical volume, i.e., 16^ lattice, /? = 6.0. The last subsection showed that for this size the 
removal of the vortices defined by the "mesonic gauge" hardly makes the string tension 
vanish. The result of 50 independent measurements is shown in figure [7| The potentials of 
this figure were fitted by the function 

V{r) 0? = a n — — , r = na . (59) 



In either case we find 



vortex only 
vortex removed 
full 



n 



aa^ = 0.061 a = 0.096 
aa^ = 0.0 a = 0.43 , 

aa^ = 0.051 a = 0.25 . 



A small Coulomb part survives the projection onto vortices. One also observes that the po- 
tential obtained from configurations from which the Laplacian vortices have been removed 
is perfectly fitted by a Coulomb law. The string tension from vortex configurations is a bit 
higher than the full string tension. This is probably due to the small physical size of the 
lattice: Coulomb contributions are dominant and are represented by the vortex matter as 
string tension to some extent. 

It turns out [TJ 121] in the case of the SU{2) gauge group that preconditioning with the 
Laplace center gauge and subsequent implementation of the "mesonic" gauge by iteration 
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Figure 8: The continuum limit of the vortex- and the monopole density for the 
"mesonic" gauge and the Laplacian gauge (left panel). The vortex- and the Z3 monopole 
density in physical units for the "mesonic gauge" (right panel). 

over-relaxation strongly reduces the influence of the Gribov copies. Center projection of 
these ensembles yields high-quality vortex matter the properties of which nicely extrapo- 
lates to the continuum limit. We have checked that repeating this approach for the SU{3) 
case produces vortex matter where the corresponding string tension again only reaches 
62% of the full string tension. 

5 Properties of SU(3) vortex matter 

In this section, we will investigate which definition of the vortex matter (MCG or LCG) 
admits an interpretation of the vortices as physical degrees of freedom. To this aim, the 
continuum limit a — of the planar vortex density Pvor (i-e., the "microscopic" vortex 
density of subsection 14. 2p and the density pmon of monopoles are investigated. 

The vortex density 

Pvor can be extracted from the ensembles {Z^(x)} as follows: we say 
that a Z3 plaquette, 

v{p) =l[Zi, / = {x,/i}, v{p)eZ,, (60) 

carries a non-trivial center flux if 

v{p) = z , z = expji^ (61) 
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and ip = —1 OT ifi = 1. li Vj3 denote the probability that a particular plaquette of the 
lattice carries a non-trivial center charge, the vortex density is defined by 

Pvor a^P) = Vp . (62) 

As outlined in subsection 13.31 the monopole world lines are associated with dual links 
whose corresponding monopole current m^(a;) (jH^ is non-zero. Alternatively, 3d hyper- 
cubes from which non-trivial center fiux emerges are said to contain a center monopole. 
Given that Qp is the probability that a particular hypercube contains a center monopole, 
the monopole density is obtained from 

Pmona\(3) = Qp. (63) 

Both quantities characterize the vortex matter. In order to interpret the vortices as sensible 
degrees of freedom in the continuum limit a ^ 0, the quantities must obey 

lim Pvor = constant , lim pmon = constant . (64) 

/3— »oo /3— >oo 

Our numerical findings are summarized in figure |H1 Figure ISJ left panel shows the densities 
Pvor and pmon in units of the lattice spacing as a function of a \/a. Simulations were 
performed for [3 = 5.6, 5.7, 5.8, 5.9. The corresponding size of the lattice spacing can be 
found in table |21 

Let us firstly focus on vortex matter obtained after implementing the "mesonic" gauge 
condition (|^. There, the data are perfectly fitted by 

Pvor a\P) ^ 0.45 (ay/^y , (65) 

Pmo„a'(/5) ^ 0.35 (ay/^y . (66) 

These findings suggest that the planar areal vortex density as well as the monopole 
density properly extrapolate to the continuum limit (see figure |Hl right panel). The same 
result for the vortex density was found for the case of a SU{2) gauge group [El E]- In a 
naive random vortex model, one expects that the string tension is given by (see 1)511) ) 

ae/ ~ ^ Pvor ~ 0.67 a . (67) 

We point out that the naive string tension 0"^/ roughly agrees with the string tension 
measured from center projected configurations (see subsection 14. 3p . i.e., 

az3 ~ 0.62 a . 

This indicates that the intersection points of the "mesonic" center vortex clusters do not 
contain significant long-range correlations. Using y/a = 440 MeV as a reference scale, one 
finds for the SU(3) case 

Pvor ^ 2.2 , p^on ^ 3.7 -\ , ^ ^ 1.16 . (68) 
fm fm 
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Figure 9: (line) density of monopoles on the shortest connected hne joining all monopole 
sites: lattice simulation (full symbols) and after randomly re-distributing the monopoles 
(open symbols). 

The latter quantity might be of interest for the construction of Z3 random vortex models. 

As already noticed for a SU{2) gauge group [2, the situation drastically changes for the 
Laplacian center gauge. The vortex and monopole densities (times the canonical powers 
of the lattice spacing) scale linearly with the lattice spacing. Two fits represent the areal 
vortex density almost equally well: 

/ \ 0.66 

Pvor a (/9) ~ 0.1 + 0.52 a .Ja , pvor a (/?) ^ 0.56 (^a , (69) 

where the latter fit function is slightly preferable. In addition, the monopole density is well 
represented by the linear function 

Pmon a' (/5) ^ 0.38 a v^. (70) 

Both quantities, i.e., pvor and pmon, diverge in the continuum limit a — 0. However, one 
finds that the ratio 

^ ^ 0.9 (LCG) (71) 

Pvor 

is roughly independent of the lattice spacing a, as is the case for the "mesonic" gauge. 

It is interesting that the monopole density (of the Laplacian gauge) diverges in a some- 
what controlled way: a situation where the monopoles lie dense on 2d hypersurfaces of 
the 4d space-time would correspond to the observed scaling with the lattice spacing. In 
order to get a rough idea of how the the monopoles are organized within space, a closed 
loop which joins all monopole sites is calculated. This loop is obtained by a simulated 
annealing procedure which minimizes the length of the loop. Thereby, the Euclidean norm 
(where in addition the toroidal topology is taken into account) is used as a measure for 
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lengths. From a mathematical point of view, finding the global minimum is a "traveling 
salesman" problem in three dimensions, and is beyond the reach of numerical calculation. 
However, the "simulated annealing" algorithm generically generates paths the length of 
which is within a few percents of the minimal length. This suffices for our purposes here. 
Given a finite set of points, it is difficult to tell whether the points are falling on top of a 
"smooth" curve. In order to gain first insight, we have calculated a connecting loop for a 
given set of Z3 monopoles with the help of simulated annealing. Dividing the number of 
monopoles by the length of this loop gives an effective line density. Finally, this quantity 
is averaged over several lattice configurations. In order get a clue about the significance 
of this average line density, we have randomly re-distributed the monopoles of a partic- 
ular configuration and we have re-calculated the effective line density. If the monopoles 
produced by the lattice simulation tend to fall on top of a smooth line, the average line 
density must be significantly larger than in the case of the random distribution of the same 
amount of monopoles. This is indeed the case, as is shown by figure El 

6 Conclusions 

MCG vortices of the gauge group SU (2) have been realized as sensible degrees of freedom 
in the continuum limit [Hj, and they are closely related to confinement [Sj. In the present 
paper, the MCG vortex matter has been investigated for the important case of SU{3), 
using large-scale numerical simulations. Focal points were the questions: To which extent 
are the SU{3) vortices relevant for confinement? Are the SU{3) vortices meaningful in the 
continuum limit? 

In a first step, we verified that the phase of large Wilson loops is represented by the 
center fiux going through the Wilson loop, we confirmed the conjecture in [211 that the Z3 
fiuxes going through a planar area are strongly correlated at length scales smaller than the 
hadronic one. As a byproduct, the fiux correlation length (see (j53|) for a proper definition) 

Li/2 ^ 0.8 fm 

was seen to be in rough agreement with scaling. 

In a second step, the static quark potential produced by the MCG vortices was addressed. 
In contrast to the case of SU{2), the string tension from the vortex projected ensembles 
turned out to be 62% of the full string tension. This finding is rather independent of the 
lattice size and the value of the lattice spacing. On the other hand, removing the Z3 vortex 
degrees of freedom "by hand" from the full lattice configurations (see flH^). always results 
in a vanishing string tension (if the physical lattice volume is large enough). This implies 
that there is still a certain relation between the MCG vortices and confinement. 

On the other hand, it is known fI^ that vortices which are defined from the Laplacian 
center gauge (LCG) reproduce the string tension to full extent. Here, we checked that 
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"preconditioning" the lattice configurations with LCG and subsequent MCG fixing does 
not produce vortex matter which yields significantly more than 62% string tension. 

The question arose as to which definition (MCG or LCG) of the vortices produces vortex 
structures which are sensible in the continuum limit. Here, the planar vortex density Pvor 
(the density of points where the vortices intersect a 2D planar hypersurface) as well as 
the (volume) density pmon of monopoles were studied, we found that in the case of the 
MCG vortices both quantities properly extrapolate to the continuum: 

^ ^ 0.45 , ^ ^ 0.35 , ^ ^ 1.16 . 

pvor 

In contrast, both quantities diverge in the continuum limit for the case of LCG vortices. 
Surprisingly, the LCG vortex matter satisfies simple scaling laws: 

Pmona=^(/3) ^ 0.38 a v^, ^ ^ 0.9 (LCG) . 

Pvor 

The investigations of sets of Z3 monopoles residing within the spatial hypercube indicated 
that the LCG monopoles tend to fall on top of a "smooth" Id curve which is embedded in 
this hypercube. 

In summary, only the maximal center gauge allows for a direct interpretation of the vortices 
in the continuum limit of pure 5'?7(3) gauge theory. There is no string tension without 
MCG vortices, but MCG vortices only support 62% of the full value. 
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